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ABSTRACT 

We present the results of the analysis of two Chandra observations of Circinus X-1 
performed in 2007, for a total exposure time of ~50 ks. The source was observed with 
the High Resolution Camera during a long X-ray low- flux state of the source. Cir X-1 
is an accreting neutron-star binary system that exhibits ultra-relativistic arcsec-scale 
radio jets and an extended arcmin-scale radio nebula. Furthermore, a recent paper 
has shown an X-ray excess on arcmin-scale prominent on the side of the receding 
radio jet. In our images we clearly detect X-ray structures both on the side of the 
receding and the approaching radio jet. The X-ray emission is consistent with being 
from synchrotron origin. Our detection is consistent with neutron-star binaries being 
as efficient as black-hole binaries in producing X-ray outflows, despite their shallower 
gravitational potential. 
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1 INTRODUCTION 

Cir X-1 is an exotic X-ray binary discovered by Margon et 
al. (1971) which shows flares with a period of 16.55 days, ob- 
served first in the X-ray band (Kaluzienski et al. 1976) and 
then in the infrared (Glass 1978; Glass 1994), radio (Haynes 
et al. 1978) and optical bands (Moneti 1992): this fact is in- 
terpreted as enhanced accretion close to the periastron pas- 
sage of a highly eccentric binary orbit (e ^ 0.8, Murdin et 
al. 1980, Nicolson, Glass & Feast 1980). Beyond variability 
at the 16.55 day orbital period the source shows dramatic 
evolution of its X-ray luminosity, spectra and timing prop- 
erties on timescales from milliseconds to decades (Shirey et 
al. 1999, Parkinson et al. 2003). The evidence that the sys- 
tem harbours a neutron star comes from the detection of 
type-I X-ray bursts (Tennant et al. 1986a,b), confirmed by 
the recent detection of twin kHz QPOs in the X-ray power 
density spectra (Boutloukos et al. 2006). Based on the prop- 
erties of the type-I X-ray bursts, Jonker & Nelemans (2004) 
estimated a distance to the source of 7.8-10.5 kpc; according 
to a different measure of the galactic absorption, laria et al. 



E-mail: p.soleri@uva.nl 



(2005) derived a lower distance of 4.1 kpc (but see Jonker, 
Nelemans, Bassa 2007). 

Cir X-1 is the most radio- loud neutron star X-ray binary 
(Whelan et al. 1977, Haynes et al. 1978; Cyg X-3 and SS433 
are brighter in radio but their nature is unclear), showing 
extended structures both at arcmin and arcsec scale. The 
arcmin-scale structure has been extensively studied by Stew- 
art et al (1993) and by Tudose et al. (2006): the source shows 
two radio jets (south-east and north-west direction) embed- 
ded in a large scale, diffuse radio nebula. There is general 
agreement that this nebula is the result of the radio lobe 
inflated by the jets over several hundred-thousand years. 
Arcmin-scale jets are curved which might be due to an in- 
teraction with the interstellar medium (ISM) or a precessing 
jet. The arcsec-scale structure has been investigated for the 
first time by Fender et al. (1998) and later studied by Fender 
et al. (2004), that reported a one-sided highly- variable jet 
(for a detailed study of the arcsec- variability evolution over 
10 years see Tudose et al. 2008), the most relativistic one de- 
tected so far within our Galaxy, with a bulk Lorentz factor 
^ 10. 

The X-ray light curve of Cir X-1 is highly variable and 
characterized by bright flares and periods of very low X- 



2 P. Soleri et al. 







■ 





Time (days from 2005 March 31st) 



Figure 1. RXTE/ASM daily-average light curve from 2005 
March 31 to 2007 November 29. The vertical dashed lines corre- 
spond to two Chandra/ imC-l observations, the vertical dashed- 
dotted line corresponds to the HETGS observation analysed in 
Heinz et al. (2007). 

ray flux (Parkinson et al. 2003). On 2005 June 2nd Cir X- 
1 was observed with the High Energy Transmission Grat- 
ings (HETGS) onboard Chandra for 50 ks during one of its 
long-term X-ray low-flux states (an example of X-ray light 
curve with long periods of low flux is reported in Figure [1]) . 
Analysing the HETGS observation, Heinz et al. (2007) found 
evidence for an arcmin-scale X-ray structure prominent on 
the side of the receding radio jet. While for black hole candi- 
dates extended X-ray jets have already been detected with 
Chandra and XMM-Newton in a number of sources (e.g. 
XTE J1550-564, Corbel et al. 2002; 4U 1755-33, Angelini 
& White 2003; H1743-322, Corbel et al. 2005), this is the 
flrst detection of X-ray structures in a secure neutron star 
system, showing that neutron stars can be as efficients as 
black holes in producing X-ray outflows. 
Here we present recent Chandra observations of Cir X-1 
where we clearly detect an extended X-ray structure both on 
the side of the approaching and the receding radio jet, con- 
flrming and extending the detection of Heinz et al. (2007). 



2 OBSERVATION AND DATA ANALYSIS 

Cir X-1 was observed with the High Resolution Camera 
(HRC; Zombeck et al. 1995, Murray et al. 1998) onboard 
the Chandra X-ray Observatory on 2007 April 21st (43 ks) 
and on 2007 May 16th (7 ks), during an exceptionally long 
interval of very low X-ray flux that occurred from ^ 2007 
February until ^ 2007 August. Figure [T] shows the X-ray 
light curve of Cir X-1 taken with the Afl Sky Monitor (ASM) 
instrument onboard the Rossi X-ray Timing Explorer satel- 
lite (RXTE) for the period 2005 March - 2007 November, 
where several episodes of low-flux can be identifled. Two ver- 
tical dashed lines correspond to our two HRC-I observations 
(hereafter observation A and observation B, respectively). 

2.1 X-ray jets 

Images from observations A and B are reported in Figure O 
The images have been rebinned and smoothed using a gaus- 
sian kernel of three pixels in radius and countour lines have 
been applied. A visual inspection of the images clearly re- 
veals the presence of an extended X-ray structure around the 



source up to '-^arcmin scale, visible in the South-East and 
in the North- West quadrant (i.e. aligned with the South- 
East - North- West direction), consistent both with the re- 
sults recently reported by Heinz et al. (2007) and with the 
ahgnment of the arcmin-scale jet observed by Stewart et al. 
(1993) and by Tudose et al. (2006); this issue will be dis- 
cussed in ^2.21 

The extended emission is detected both in obs. A and B 
and the two images (note the different exposure in the two 
observations), at a visual inspection, are consistent: the X- 
ray emission around the central source, in both of them, 
elongates along the same axis and furthermore, it presents 
similar structures in the North- West quadrant and a simi- 
lar spike-like structure in the northern quadrant (although 
there are differences specially in the South-East quadrant, 
due to different sensitivity limits) . 

Since the two images are consistent, we can merge them to 
have an image with a higher equivalent exposure, using the 
standard Ciao 4-0 analysis tools (Fruscione et al. 2006). We 
adaptively smoothed the image obtained from the merging 
with an average signiflcance of 8.3a per smoothing length 
(the minimum smoothing length is 6 pixels, the maximum 
one corresponds to the image size) and we applied contour 
levels. The resulting image is plotted in Figure [2] (panel (c)), 
where a diffuse X-ray emission elongated along the South- 
East - North- West direction is evident, extending up to ^1 
arcmin from the point source. Other structures might look 
real from an analysis of this image (e.g a circular excess in 
the South- West quadrant) but they will not be discussed 
since an inspection of the residuals obtained after fltting the 
point spread function (PSF) suggests that they are noise 
features (see next part of this paragraph). 
Figure [3] shows two proflle cuts (for obs. A), one extracted 
on a region aligned with the X-ray excess and one along the 
perpendicular axis, showing an excess along the South-East 
- North- West axis and supporting the evidence that there is 
an X-ray extended structure aligned with a privileged direc- 
tion. 

From Figure [2] (panels (a), (b) and (c)) and Figure [3] the 
presence of an extended X-ray structure around the central 
source is evident but to make our detection more robust, an 
analysis of the PSF is needed. We simulated a monochro- 
matic PSF at 1 keV (from Heinz et al. 2007 we expect the 
X-ray excess emission to peak at this energy) using ChaRT 
(Carter et al. 2003), considering the same number of counts 
as detected from the source and we projected it in the de- 
tector plane using the Marx ray tracing (Wise 1997). Since 
the PSF wings depend on the position on the detector, for 
the PSF analysis we used only the obs. A: even if the im- 
age from the merging has a higher equivalent exposure, the 
PSF extracted from that image might contain artefacts (the 
source position on the detector is not the same for obs. A 
and B) that we want to avoid. For the PSF analysis, the 
Sherpa 3.4 tools (Freeman et al. 2001) have been used. 
In Figure [2] (panel (d)) we show the image of the residuals 
obtained after fltting the source image with a gaussian func- 
tion and a constant using the PSF image as the convolution 
kernel. The presence of an extended X-ray structure aligned 
with the South-East - North- West axis is evident, at a dis- 
tance from the central source between ~ 25" and ^ 50" (the 
lower limit is inferred by inspecting both the residuals image 
and the radial proflle image in Figure [3]). The X-ray excess 
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Figure 2. (a) top left: HRC-I image of obs. A. (b) top right: HRC-I image of obs. B. (c) bottom left: image of the merger of obs. A 
and B, adaptively smoothed, (d) bottom right: image of the residuals obtained after fitting the source image (obs. A) with a gaussian 
function and a constant using the PSF image as the convolution kernel. In all the panels we marked the source location with a cross. 



on the side of the receding radio-jet (North- West quadrant, 
hereafter "receding X-ray jet") lies at position angle (PA, 
measured counterclockwise from due North from the point 
source) intervals 286°-295° and 307°-324°. A knot appears 
also at PA ~ 10° and its identification will be discussed in 
the next sub-section. On the side of the approaching radio 
jet (South - East quadrant) we see one main X-ray blob 
in a PA interval 88°-152° (elongated along a North-East - 
South- West direction, hereafter "approaching X-ray jet"), 
at a distance from the point source between 28" and 38". 
Other minor blobs appear in the residuals image: all of them 
are located in the South - East quadrant and are consistent 
with being knots of the approaching X-ray jet. 

2.2 Comparison with previously detected X-ray 
extended structure and the radio jet 

The radio nebula and the relativistic radio-jets of Cir X-1 
have been investigated during multiple epochs of radio ob- 
servations, both on arcmin and arcsec scale (Stewart et al. 
1993, Fender et al. 1998, Fender et al. 2004, Tudose et al. 
2006, Tudose et al. 2008). From those papers, we estimate 
the corresponding PAs for the arcmin-scale jets: 303° =b 9° 
and 310° =b 15° for the jets in Stewart et al. (1993) and 
Tudose et al. (2006) respectively and 320° =b 4° for the ul- 
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Figure 3. Radial surface-brightness profile extracted, for obs. 
A, across the jet (circles) and across a region perpendicular to 
the jet (stars). The extraction regions are shown in the inset and 
correspond to the following PAs (see definition in the text): X-ray 
excess 106°-144° and 286°-324°; axis perpendicular to the X-ray 
excess 16°-54° and 196°-234°. 

trarelativistic jets in Fender et al. (2004). All these angles 
are consistent with the PAs of the X-ray jets in our HRC-I 
observations and fall between the two X-ray filaments de- 
tected by Heinz et al. (2007) in the HETGS image. 
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Figure 4. Radio- X-ray overlay. Grey scales: X-ray residuals im- 
age (Fig. [2] (d)); green contours: adaptively smoothed image 
of obs. A &; B merged (Fig. [2] (c)); red contours: adaptively 
smoothed, normalized, PSF subtracted image (from Heinz et al. 
2007); dark blue contours: 1.4 GHz surface brightness (adapted 
from Tudose et al. 2006, levels increase by \/2 between contours, 
outermost contour: 11.2 mJy/beam); light blue lines: estimated 
allowed range of PAs from high-resolution radio observations of 
approaching radio-jet (Fender et al. 2004). The crossing of these 
two lines identify the source position. 

Figure |4] shows an overlay of the residuals image (the same 
as Figure [2] panel (d), grey scale), the X-ray contours from 
Figure [2] panel (c) and the jet emission from Figure 1 of 
Heinz et al. (2007). Also shown are the radio contours from 
Figure 3 of Tudose et al. (2006) and the limits on the PA for 
the arcsec-scale jet (Fender et al. 2004). The X-ray jets de- 
tected in the HRC-I images are broadly consistent with the 
X-ray excess of Heinz et al. (2007). Besides the X-ray knot 
that we detected at PA 10°, the consistency between the 
X-ray and the radio jets is clear on both sides (approaching 
and receding) : what we see in the HRC-I images appears as 
the X-ray counterpart of the radio jets from Cir X-1. 



3 DISCUSSION AND CONCLUSIONS 

The analysis of our two HRC-I observations of Cir X-1 
clearly showed X-ray jets both on the side of the receding 
and the approaching radio jet, consistent with being the X- 
ray counterpart of the arcmin-scale radio jets (Tudose et al. 
2006). Heinz et al. (2007) proposed two possible alternative 
explanations for the origin of the X-ray excess: synchrotron 
emission and thermal bremsstrahlung. We now investigate 
whether these two mechanisms can still explain the X-ray 
jets as observed by HRC-I. 

Synchrotron emission: HRC-I does not have the en- 
ergy resolution necessary to allow for spectral fitting, in- 
stead we calculate the jet flux with webPIMM^ usmg, as 
input spectrum, the absorbed power law fitted in Heinz et al. 
(2007) (F = S.Oli s, = 5.9 x 10^^ cm"^) and the average 
count rate extracted, for obs. A, in the region correspond- 
ing to the main knot of the approaching jet, using Ciao 4-0. 

^ http:/ /heasarc. nasa.gov/Tools/w3pimms. html 



This gives an un-absorbed flux F2-iokeV = 7.4 x 10~^^ erg 
cm~^ s~^, corresponding to an X-ray luminosity L2-iokeV — 
5.4 X 10^^ erg s~^D^_q (throughout the paper we adopt a 
distance from the source D = 7.8 kpc, the same adopted 
in Heinz et al. 2007, to make comparisons easier). The used 
photon index F is broadly consistent with the emission be- 
ing of synchrotron origin {J^syn ^ 1*5)5 Q^ssuming a source 
volume V, a specific luminosity L^^ and a spectrum of the 
form (X v^, (a = 1 — F), we can follow Longair (1994) 
and Fender (2006) to estimate the minimum energy associ- 
ated with the source. The morphology of the jet suggests 
that we are observing the surface of a conical volume where 
the emission takes place. Considering the jet length to be 
roughly I jet ~ 1-6 pc D7.s/sini (where i is the angle with 
the line of sight) and a half-opening angle [3 ~ 19.3° the vol- 
ume of the emitting cone is Vjet = 9.7 x lO^^cm^i^f.g/smi, 
assuming equal emitting volumes both on the receding and 
on the approaching side. Under these assumptions the mini- 
mum jet energy is Emin ^ 6 x 10^^ erg. Again following Lon- 
gair (1994) and Fender (2006), we can calculate the magnetic 
field Bmin associated to Emin, the Lorentz factor 7 of the 
energetic electrons emitting by synchrotron and their gyro- 
radius rg: Bmin = 8.2/iG, 7 = 2.3 x 10^ and rg = 1.6 x 10"^ 
pc. The energy of these electrons is = 1.2 X 10^^ eV and 
their energy loss rate is R ^ 2300 eV/s: considering that we 
have no evidence for re-acceleration taking place, the life- 
time tsyn of these electrons is tsyn ~ 1600 yr. Following 
Tudose et al. (2006), we assume that the jets are injected 
mainly during the flare states (duty cycle ^^6%): the result- 
ing minimum jet power is Wsyn ^ 2 x 10^^ erg/s. 

Thermal bremsstrahlung: in this case the X-rays 
would originate from the shock driven into the ISM by the 
propagation of the jets. Here we use a model developed 
for extragalactic jet sources by Castor et al. (1975), Kaiser 
& Alexander (1997) and Heinz et al. (1998). We assume 
the temperature of the thermal gas to be kBTshock = 2.2 
keV (kBTshock = 2.21^-; keV and Nh = 5.4 x 10^^ cm^^ 
are obtained by Heinz et al. 2007 fitting the jet spectrum 
with an absorbed thermal model), the electron density to 
be He ~ 10 cm~^ (as used in Heinz et al. 2007; since the 
temperature Tghock = 25.5 MK we assume an ionization 
fraction x=l: all the hydrogen is ionized) and the length of 
the shock region (considering both receding and approach- 
ing side) Lshock ~ 2.27 pc. Such length and the used den- 
sity give an emitting gas mass of 0.02-0.4 MqD^^q /sini: 
uncertainties are due to different possible measures of the 
thickness of the shock region. Following Kaiser & Alexan- 
der (1999) and balancing the interior pressure exerted by the 
jets and the ram pressure of the shocked ISM, the jet lifetime 
is tth = f ~ I'^'OO yr (vsp is the velocity of the shock- 
compressed particles, obtained from the shock temperature 
ksTshock = 2.2 keV) and the jet power is Wth ^ 6 x 10^'^ 
erg/s. 

Discussion: In the synchrotron case the jet lifetime 
tsyn is smaller than the time expected for the jet to inflate 
the large-scale radio lobe (lO'^ — 10^ yr, Tudose et al. 2006) 
and this suggests that the X-ray emission could come from 
the jet itself rather than from the inflated radio nebula. Fur- 
thermore, the gyro-radius of the electrons responsible for the 
X-rays is smaller than the size of the jet (rg < I jet) and this 
implies that those electrons can be confined in the jet re- 
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gion. Our jet power value is consistent with the estimate of 
Tudose et al. (2006) (Wjet.Tudose ^ 10^^ erg/s, which how- 
ever could be up to two order of magnitude higher) and one 
order of magnitude smaller than the jet power calculated 
by Heinz et al. (2007) {Wsyn, Heinz > 5 x 10^^ erg/s). The 
magnetic field associated to the emitting electrons is again 
consistent with the value reported in Tudose et al. (2006) 
{Bmin,Tudose = 6.3/xG), making Wsyn a robust estimate of 
the jet power, sufficient to inffate the radio nebula. Wsyn 
is > 0.1% of the Eddington luminosity {LEdd = 1.8 x 10^^ 
erg/s) for a IAMq neutron star: Cir X-1 is only slightly 
super-Eddington and this would imply a jet-production ef- 
ficiency r] > 0.01%, consistent with what obtained by Heinz 
et al. 2007 (ry > 0.5 %). 

In the thermal bremsstrahlung case we calculated a jet 
power Wth ~ 6 X 10^^ erg/s, two orders of magnitude 
larger than Wsyn and possibly not consistent with the val- 
ues estimated by Heinz et al. (2007) {Wth, Heinz = 5 x 10^^ 
erg/s). Furthermore, Wth is a significant fraction (^ 33%) 
of the Eddington luminosity {Lsdd — 1.8 x 10^^ erg/s) 
for a IAMq neutron star and this would be an extremely 
high jet-production efficiency 77 ^ 3%, even for an accreting 
black hole. Taking an Eddington-limit mass accretion rate 
M = 10^^ g s-^ and even considering that all the accreted 
mass is ejected in an outfiow, the time required to infiate 
the jets would be ~2.7 Myr, much bigger than ttn ~ 1700 
yr, suggesting that the emission can not come from the jet 
itself. Therefore thermal bremsstrahlung appears as an un- 
likely emission mechanism (especially compared to the syn- 
chrotron case) for the X-ray jets. 
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